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METHOD FOR THE LOCALIZATION OF ONE OR MORE TRANSMITTERS 

The invention relates to a method for the localization of one or 
more transmitters. 

5 It can be applied essentially to the localization of transmitters on 

the ground from a moving vehicle without a priori knowledge concerning the 
signals sent. 

The technical field in particular is that of the passive localization of 
transmitters. 

10 Figure 1 illustrates an airborne localization. The transmitter is at 

the position (xo,yo,zo). The carrier at the instant t k is at the position (x^y^ f z k ) 
and perceives the transmitter at the angles of incidence 
WfcXo^o.Zo), A{t kt x 0i y 0$ z 0 )). The angles ^f,x 0 ,y 0 ,zo) and A(t,x 0i y 0 ,zo) evolve in 
time and depend on the position of the transmitter as well as the trajectory of 

15 the carrier. 

The angles 6(f,x 0 ,yaZo) and A(t,x 0 ,yo,Zo) are determined relative to 
a network of N antennas that can be fixed beneath the carrier as shown in 
figure 2. 

There are several existing techniques to determine the position 
20 (x m , y m , z m ) of a transmitter. These techniques of localization differ 
especially in the parameters which are estimated instantaneously at the level 
of the network of sensors. They can be classified as follows. 

Use of goniometry 

25 These techniques are known and used in the prior art. In most 

cases, they are based on a 1D azimuthal goniometry. The azimuths 
0km=6(t k ,x mi y mt z m ) associated with the m th transmitter are measured for 
different instants t k . In using the position (x k ,y k ,z k ) of the carrier at the 
corresponding instant k, a position (x m ^,y m ^,z m ^) of the transmitter m is 

30 estimated by a ground intersection. The position (x k ,y k ,z k ) of the carrier is 
given by a GPS unit, its orientation is obtained by a compass in the case of a 
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ground carrier and by an inertial navigation system in the case of an aircraft. 
From all the positions (x mk ,y mkf z mk ) t the method extracts data with which it is 
possible to determine the M dominant positions (x m ,y m ,z m ) of the incident 
transmitters. The localization is obtained by triangulation or by ground 
5 intersection (2D goniometry). The drawback of triangulation techniques is 
that they require major movement. Furthermore, goniometry techniques 
require the use of a network of non-ambiguous sensors. This has the 
drawback of necessitating a calibration table and restricting the size of the 
sensor network and therefore giving incidence values of limited precision. 

10 

Use of the phase difference between two remote sensors 

The inter-sensor phase difference A(p(t k ,xo,yo,zo) depends on the 
positions of the two sensors as well as the incidence 

15 fl(fo,xo,yo,zo), A(t k ,x 0 ,yo,zo)) of the transmitter. This phase, which depends on 
time, is directly related to the position (x 0l yo,zo) of the transmitter. 
Consequently, studying the function of time A<p(t,Xo,yo,Zo) makes it possible to 
deduce the position (xo,yo,zo) of the transmitter. In this group of applications, 
the two sensors are distant in order to augment the precision of 

20 measurement of the phase. This has the drawback of causing variations in 
the phase difference A<p(t,Xo,yo,Zo) as a function of time over more than 2n 
and the technique then necessitates a step enabling the phase to be made to 
vary by more than 2n. Furthermore, in this technique, the phase is measured 
by carrying out a direct intercorrelation between two sensors, and cannot be 

25 used to deal with the multiple-sensor case. 

Use of the measurement of the carrier frequency of the transmitter 
These techniques make use of the fact that the estimated carrier frequency is 
the sum of the carrier frequency of the transmitter and the Doppler shift due 
to the speed of movement of the carrier. The Doppler shift has the 

30 advantage of depending on the position (x 0 ,yo,Zo) of the transmitter and of 
also being a function of the time Af(f,x 0 ,yo,z 0 ). Consequently, studying the 
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function of the time Af(f,x 0 ,yo,zo) makes it possible to deduce the position 
(x 0j yo,2o) of the transmitter therefrom. However, the measurement of this 
Doppler shift has the drawback of necessitating transmitters with particular 
waveforms. This measurement of frequency can be done by cyclical 
5 techniques in which it is assumed that the signal sent is non-circular. 
Use of propagation times 

These techniques use the differences in propagation time between 
antennas (TDOA or time difference of arrival) which are directly related to the 
respective differences between the transmitter and the different antennas 

10 and therefore to the position (x 0l yo,z 0 ) of the transmitter. The use of at least 
three antennas that are sufficiently spaced out enables the position (x^yo.zo) 
of the transmitter to be deduced by hyperbolic localization. The drawback of 
these techniques is that they cannot be implemented in a single-carrier 
context owing to the considerable spacing required between antennas. 

15 Furthermore, in these techniques, the time difference is measured by the 
direct performance of an inter-correlation between two sensors. This 
approach cannot be used to deal with the case involving multiple 
transmitters. 

The method of the invention relies especially on a novel approach 
20 for the direct estimation of the positions (x m , y m , z m ) of each of the 
transmitters from a parametric analysis of the multipath signal at various 
instants t k on a duration Dt. The parametric analysis especially has the 
additional function of separating the different transmitters at each point in 
time t k . Then, the parameters of a same transmitter coming from the different 
25 points in time t k are associated to finally localize each of the transmitters. 

The invention relates to a method of localization of one or more 
sources, said source or sources being in motion relative to a network of 
sensors, the method comprising a step of separation of the sources in order 
to identify the direction vectors associated with the response of the sensors 
30 to a source having a given angle of incidence. The method is characterized 
in that it comprises at least the following steps : 
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• associating the direction vectors a^.-.a^ of the m th transmitter 
obtained respectively at the instants t r ..t K , 

• localizing the m th transmitter from the associated vectors a^.-.a^. 

5 The method of the invention has especially the following advantages : 

- it enables the localization of the position of the transmitter in (x, y, z) 
coordinates and, in addition, that of its speed vector, 

- it can be applied when there are one or more incident transmitters 
present, 

10 - its implementation does not require any particular knowledge of the signal 
sent, 

- it enables the use of a network of ambiguous sensors, i.e. several 
incidence values are associated with the same response of the network 
which have the advantage of being great and therefore of being more 

15 robust with respect to inter-antenna coupling phenomena or, more 
generally, antenna-network modeling errors, 

- it can be implemented in networks calibrated in (6, A), 

- it can be implemented in networks with amplitude-diversity antennas such 
as co-localized antennas : networks with dipoles having the same phase 

20 center and different orientations. 

Other features and advantages of the object of the present 
invention shall appear more clearly from the following description given by 
way of an illustration that in no way restricts the scope of the invention 
and from the appended figures, of which : 
25 Figure 1 shows the principle of the localization by an aircraft of a 

transmitter positioned on the ground, 

Figure 2 shows the relationship between a network of antennas 
and the value of incidence of a transmitter, 

Figure 3 is a general drawing explaining the working of the method 
30 according to the invention, 
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Figures 4, 5 and 6 are examples of implementation of the method 
according to the invention. 

In order to provide for a clearer understanding of the object of the 
present invention, the following description is given by way of an 
5 illustration that in no way restricts the scope of the invention for the 
localizing of several transmitters positioned on the ground using a 
network of sensors installed in an aircraft in motion. Such a system is 
described for example in figure 1. The aircraft is equipped with a 
processor adapted to implementing the steps of the method according to 
10 the invention. 

The method can also be implemented in the context of vehicles 
moving on the ground. 

Figure 3 is a representation, in a signal time/amplitude graph, of 
the signal x(t) formed by a combination of the signals of the transmitters 
15 at different instants ti, t 2 , ...tK. This figure summarizes the different steps 
implemented, namely the separation of the transmitters SE and the 
parametric estimation EP, the association of the parameters of each 
transmitter and the localization of a transmitter. 

In the presence of M transmitters, at the instant t at output from 
20 the N sensors of the network, the method has the vector x{t) representing 
the mixture of signals from the M transmitters. Around the instant t k , the 
vector x(t+t k ) sized Mx1 , representing the mixture of the signals from the 
M transmitters, is expressed as follows : 

x(f+fc)=i; a(0* m ,A to )s m (/H^ pour |f|<Af/2 

m=l 

(D 

25 where b(t) is the noise vector assumed to be Gaussian, a(9,A) is the 
response of the network of sensors to a source having an incidence (9,A), 
A/FlafG/rt, A M )... a(9* M , A**,)], s{t)=[s^t)...s^t)] J t Q km =0{t k ,x mi y mi z m ) and 
Afrm^ffcX/n.ym.Zm). In this model, the mix matrix A* depends on the instant t k 
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of observation. 

The direction vector with incidence corresponding to the m th transmitter at the 
instant t k 



is a known function of t k and of the position of the transmitter (x mi y mt z m y 



The method according to the invention comprises at least the 



following steps: 

1 . estimating one or more parameters associated with the position of the 
source, for example the direction vectors, the incidence values, the 
position etc. and separating the M transmitters for the different instants t|<: 
this consists in identifying the direction vectors having an incidence a km 
for (1<A77<M). This first step is performed, for example, by means of 
source separation techniques known to those skilled in the art. 

2. associating the parameters estimated for the m th transmitter, for example 
in associating the different direction vectors having incidence values , 
a<| m ak m obtained respectively at the instants t-|,...tk, 

3. localizing the m th transmitter from the associated vectors. 
Association step 



In the presence of M transmitters and after separation of sources, 



the method, at the instant t k , possesses the M signatures a km for (1<m</W). 
At the instant the source separation gives the M vectors b/ for (1</<M). 
The goal of this tracking is to determine, for the m th transmitter, the index i(m) 
which minimizes the difference between a km and b/fmj. In this case, it will be 
deduced therefrom that a k + 1t m = bj (m) . To make this association, the distance 
between two vectors u and v is defined, for example, by : 



a km = a(9* m> A/c m )=a(^,x m ,y m ,z m ) of the m th transmitter 



K2) 



d(u,v) =1- 




(3) 



Where u H is the conjugate transpose of the vector u. 
Thus, the index i(m) verifies : 
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d(a* m , bj(m)) 



min [4(a km , 6,)] 



(4) 



hence 



d(a km , bj(m)) 



min [1 




In this association, we consider a function p m associated with the m th 
transmitter : 



5 In the course of the association, for each transmitter m is obtained and for 
1<m<M, the function p m (f) is obtained. This function is aimed especially at 
eliminating the instants t k whose value $ m (t k ) appears to be far too distant 
from an interpolation of the function p m (f), i.e. the aberrant instants which 
may be associated with other transmitters are eliminated. A zone of 

10 tolerance +/-A is defined around the curve defined by the function p m (fo). This 
zone of tolerance will depend on the precision of estimation of the direction 
vectors a km . In particular, in the presence of M=1 sources, the zone will be in 
the range of A=3/V5At (where At is the elementary time of parametric 
estimation illustrated in figure 3 and B is the instantaneous band of the signal 



following : 

Step ASE - 1 : Initialization of the process at k=2. The initial number M of 
transmitters is for example determined by a test to detect the number of 
20 sources at the instant t 0 known to those skilled in the art, 

Step ASE - 2 : For 1<m<M determining the indices i(m) in applying the 
equation (4) and in using the vector a kim with 1<m</W and the vectors b, 
identified at the instant fo+i for (1</</W), 

Step ASE - 3 : For 1 <m<M performing the operation a k+1 m = bi (m)i 
25 Step ASE - 4 : Incrementing /«-/c+1 and if k<K returning to the step ASE-1 , 



Pm(fo) - d(a km , a 0m ) 



:<s) 



15 X(t)). 



The steps of this association for K instants t k are for example the 
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Step ASE - 5: Starting from the family of instants 0={ fi<...< t K }, 
eliminating the / instants f,E<P such that the coefficients p^f,) do not belong to 
the zone demarcated by the interpolation curve of the p m (tk) values and the 
zone of tolerance A. The instants t k where |pm(tk)- Mtk-i)|<A will also be 
5 eliminated. After this sorting, the new family of instants is 0={ h<...< t f } and 
we assume K=l. 

At the end of these steps, the method has determined the vectors 

ai m aKm associated with the m xh transmitter. 

Localization of a transmitter 

10 The method determines the position of the m th transmitter from 

the components of the vectors ai m up to ax™. These vectors a km have the 
particular feature of depending on the instant t k and above all on the position 
(Xm.ym.Zm) of the transmitter. In particular, for a network formed by A/=2 
sensors spaced out by a distance of in the axis of the carrier, the vector 

15 verifies a km : 

1 

Bkm " exp^y2tf j cos(0{t k ,x m9 y m ,z m )) cos( A^ , x m ,y m , z m )) 

The value 1 of the first component corresponds to the reference 
sensor. According to figure 1, the incidence (d(f/(,x m ,y mi z m ), A(t k ,x m ,y m ,z m ) ) 
can be directly computed from the position (x kf y k ,z k ) of the carrier at the 
instant t k and the position (x m ,y m ,z m ) of the transmitter. 

20 Step of transformation of the vector 

According to a first alternative embodiment, the method comprises 
a step of correction of the a kmt values, the measurement of the direction 
vectors a km is generally obtained to within an indeterminate complex factor. 
According to this first variant, the method comprises a step of reducing the 

25 phase reference of the direction vector measured in reducing the operation to 
the phase barycenter (defined to within a constant scalar coefficient that may 
be set at 1). This operation is performed, for example, by estimating the 



(6) 

= a(^,x m ,y m ,z m ) 
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correction coefficient determined by the following conversion of the values 
a km into values a' km : 




5 The correction coefficient is not totally determined by this 

expression given the Nth order indeterminacy of the complex root. A tracking 
of the phase evolution during the observation period is therefore done. 

Since the complex coefficient is defined to within a factor from 
among the N N th roots of the unit, the phase tracking consists in arbitrarily 

10 setting the first correction coefficient (in taking the root 1 for example), and 
then in determining, at each new iteration k+1, the coefficient that minimizes 
the mean phase differences between the direction vector recentered at k+1 
and the vector recentered at the instant k. 

The minimization criterion for measurement at the same frequency 

15 may be equal to : 



min £rain(mod( 

P 9 ** An* 



atg( 



2n),2n - mod( 



),2*)) 



(8) 



where the values a k +i are the direction vectors recentered with the correction 
coefficient arbitrarily determined by any one of the Nth roots of the 

20 expression. For measurement at different frequencies, it is possible to 
compare the phases of the components of the two direction vectors in 
correcting them by a power given by the ratio of these two frequencies. 

If we consider the vectors b km = a\ m , it is then possible to 
compare this measurement with the theoretical value b(t k x m ,y m z m ) for which 

25 the theoretical direction vector a(t k x m> y m z m ) is computed for a theoretical 
origin at the theoretical (geometrical) phase barycenter (geometrical locus for 
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which the theoretical sum of the phases differences gets cancelled out). This 
locus (generally) does not coincide with the phase center of the network. 

According to another alternative embodiment, the method 
comprises a step of conversion of the vector a km into a vector b km whose 
components are formed out of components of the vector a km . In particular, 
the method builds for example the vector b km sized (N-1)x1 in choosing a 
reference sensor in n=i : 



i*(*-l)/«i«(Q 

■ (' + l)/«*.(0 



»i.<tf)/«*(i) 



(9) 



= b(t k ,x m ,y m ,z m ) 



where a km {i) is the / h component of a km 

The components of b km correspond in this case to the ratios of the 
components of the vector a km and of the vector a^i). 
10 Thus in the example of the equation (6) in fixing /=1 we get : 

(10) 



exp 



jln - cos(0{t k ,x m9 y m ,z m )) cos( A(t k ,x m ,y m ,z m )) 
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= a km (2)/ a km (1) 

It being known that the direction vectors a km are estimated with a 
certain error e km such that a km - a(f/f,x m> y mj z m )+ e km t it can be deduced 
therefrom that the same is true for the converted vector b km of (9) . 
Step of maximization of a correlation criterion 

It being known that the vector a km is a function of the position 
(x m ,y m ,z m ) of the transmitter, the same is true for the vector b* m . The method 
comprises a step of maximization of a normalized vector correlation U(x,y,z) 
in the space (x,y,z) of position of a transmitter where 
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b K H v K (x,y,z)\ 



(11) 



Mx,y,2)= 



with 



lm 



h{t x ,x,y,z) 



= v K (x m ,y m ,z m ) + w K , v K (x,y,z) = 



b(^,x,j;,z) 



w 



lm 



and wk = 



£m 



The noise vector w K has the matrix of covariance R=E[w* w* h ]. 
Assuming that the matrix R is known, the criterion can be envisaged with a 
whitening technique. Under these conditions the following L/<(x,y,z) criterion 
is obtained : 



With R=E[w K w K H ] 

5 • It must be noted that the criteria of the equations (11) and (12) are equal 
when R=o 2 l, i.e. when the errors are considered to be of equal level on all the 
sensors and independent between sensors. The criterion L/^x^z) of the 
equation (12) is therefore valid for a white statistics model noise Wk, 
• The criteria L/^y.z) and l_K(x,y,z) range between 0 and 1 and verify 

10 L,<(x,y,z)= U'(x,y,z)=1 for the position (x mj y mj z m ) of the m th transmitter. This 
normalization enables the setting of an efficient threshold of localization r|. 
Thus, all the maximum values (x m ,y mi z m ) of U(x,y,z) which verify 
U(Xm,ym,z m )> r| are considered to be good localizations. The threshold may 
be set as a function of an approximate knowledge of the statistics of w K , 



b/R" 1 v K (x,y,z) 



(12) 



L K '(x,y,z)= 



(b/R- l b K lv K (x,y 9 z) H R-\ K (x,y 9 z)) 
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• The criteria U(x,y,z) and U'(x,y,z) have the advantage of being capable of 
implementing a technique of localization in the presence of a network of 
sensors calibrated in space (0,A). It being known that, at the instant fo, the 
analytic relationship linking the incidence (0(k,x,y,z), A(k,x,y,z)) of the 

5 transmitter with its position (x,y,z) is known, it is then possible, from the 
incidence (9(f/c,x,y,z), A(fo,x,y,z)), to deduce the vector a(f/c,x m ,y m ,z m )= 
a(9(f*,x,y,z), A(f*,x,y,z)) in carrying out an interpolation of the calibration table 
(relative to the calibrated antennas). It may be noted, however, that this 
method is insensitive to a phase bias (owing to the vector correlation 

10 criterion), 

• These criteria also make it possible to take account of the phase and 
amplitude of the components of a(6,A). The method can therefore be 
envisaged with networks of co-localized antennas with pattern diversity. 

It must be noted that, in an airborne context, knowledge of the altitude h of 
15 the aircraft reduces the computation of the criterion in the search space (x,y) 
in assuming z=/7. In the example of the equations (6) and (10), the vector v K 
(x,y,z) is written as follows : 

(13) 



v K (x,y,z) = 



( d ^ 
exp j In — cos(0fo , x, y, z)) cos( A(f , , x, y, z)) 



exp j In T cos(0(t K , x, y, z)) cos( A(^ , x, y, z)) 

A 



In this method, it is possible to envisage initializing the algorithm at 
K-Kq and then recursively computing the criterion L/<(x,y,z). Under these 
20 conditions, U<(x,y,z) is computed recursively as follows : 



U.,(x,y.z)= p^pL 



(14) 



where 
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a K+ i(x,y,z)= a^x.y.z) + b K+t m H b(f K+1 ,x,y,z) 
YK + i(x,y,z)= YK(x,y,z) + b(f K+ i,x,y,z) H b(f K+1 ,x,y,z) 

Pk+1 = Pk + bK+J m H bx+f m 

The coefficients a K+ i(x,y,z)= a^x.y.z) YK + i(x,y,z)= y/<(x,y,z), Pk + i=Pk 
are intermediate spectra enabling the computation of U + i(x,y,z). 

When the vectors b(fK+i,x,y,z) and b km have constant norms equal 
to p\he relationship of recurrence of the equation (14) becomes : 

Wx , y , 2)= K r <*-*f < 15 > 

J3 2 (K + 1) 2 

where 

aK + i(x t y,2)= ctK(x,y,z) + b K+ tm H b(f K+ i,x f y,2) 

5 The method has been described hitherto in assuming that the 

vectors have fixed positions. It can easily be extended to the case of moving 
targets with a speed vector (v xm ,Vy m> Vzm) for which there is a model of 
evolution. In these conditions, the incidence of the m th transmitter can be 
parametrized as follows : 

Q«m=0{t k , X m -V xm t k ,ym-V ym t k , Z m -V zm t k ) (16) 
et Akm=A( t k , X m -V xm t k ^m-Vym t k , Z^Vzm t k ) 

10 where (x m ,y m ,z m ) is the position of the transmitter at the instant t 0 and 
(Vxm.Vym.Vzm) are the components of the speed of the transmitter at the instant 
fo. In these conditions, the vector b* m of the equation (9) is parametrized by 
(x m ,y m ,z m ) and (Vx^Vy^Vz™) as follows : 

b km = b(^x m ,y m ,z m ,v xm ,v ym ,v zm ,)+ w km (17) 

Naturally, the criteria of localization L K and L K ' of the equations (11) and 
15 (12) are no longer parametrized solely by (x,y,z) but also by (v x ,v y ,v z ). The method 
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then consists in maximizing the criterion U(x,y,z, v Xl v y ,v z ) as a function of the 6 
parameters (x,y,z, v x ,v y ,v z ). 

The method can be applied to a very large number of measurements. In 
this case, the method comprises a step of reduction of the numerical complexity of 
5 computation (which depends on the number of measurements) by reducing K. The 
method provides for the following processing operations to be performed on the 
elementary measurements : 

- decimation of the instants t kf in eliminating the neighboring instants for 
which the development of the curve p m (/*) is not significant, 

10 - filtering (smoothing of the measurements which are the direction vectors) 
and sub-sampling, 

- the measurements are then merged on a defined duration (extraction by 
association of direction vector to produce a synthesis measurement). 

Recapitulation of the steps of the method 

15 The process of localization of several transmitters using K instants t k 

can be summarized by the following steps : 

Step No. 1 : Identification of the vectors a km for (1</7)<M) at the K instants f k 
in applying for example a technique of source separation, and source 
identification as described in the references [2] [3]. 
20 Step No.2 : Association of the vectors ai m up to a Km obtained at the 
respective instants t K associated with the m th transmitter for \<m<M in 
applying the steps ASE-1 to ASE-5 described here above. 
Step No. 3 : Initialization of the process at m=1 . 

Step n°4 : Conversion of the K vectors a km into vectors b Km as suggested by 
25 the equation (9). 

Step No. 5 : Computation and maximization of the criterion U(x,y,z) of the 
equation (1 1 ) to obtain the position (x m ,y m ,z m ) of the m th transmitter. 
Step No. 6 : Incrementing /m-m+l and if m<M return to the step No. 3. 
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In order to refine the estimation of the position (x m ,y mi y m ) of the transmitters, 
the steps of the method may be performed iteratively as follows: 
Step No. 7 : Identification of the vectors b, for (1 <i<M) at the instants f«+i in 
applying for example a technique of source separation and source 
5 identification as described in the references [2] [3]. 

Step No. 8 : For 1<m<M determining the indices i(m) in applying the 
equation (4) and in using the vector a km and the vectors b, for (1 <i<M). 
Step No. 9: For 1<m<M the operation a K +i m= b i(m ) is performed. 
Step No. 10 : For 1<m<M computation of the criterion U+i(x,y,z) iteratively 
10 in using the equations (14) and (15) and minimization of U + i(x,y,z) to obtain 
the position (x m ,y m ,Zm) of the m {h transmitter. 

Step No. 11 : If it is decided to go on, in order to achieve greater precision 
and reduced ambiguity, then the method returns to the step No. 7. 

15 Exemplary implementation of the method 

The simulations were made with a network of A/=2 sensors 
aligned in the axis of the carrier with dA = 3. Since d/A. = 3 a method 
performing goniometry at the instants t k , it would be completely 
ambiguous and would not subsequently enable operations of 

20 triangulation for localizing the transmitter. In figures 5, 6 and 7 
corresponding to the criteria of localization for K = 3, 7 and 16, the 
pseudo-spectrum L K (x,y) to be maximized, enabling the position of the 
transmitter in the space (x,y) to be determined, is plotted. It being 
known that if the transmitter is situated in (x 0 ,y 0 ), then L K {x 0 ,y 0 ) = 1, it is 

25 deduced therefrom that the iso-level curves L K (x,y) = 0.99 characterize 
the width of the major lobe. Observing that the precision of localization 
depends on the width of this lobe, it is deduced therefrom, according to 
figures 5, 6 and 7 that the greater the value of K , the greater will be the 
localizing precision. 

30 
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